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Abstract. The well-poscdncss of a phase-field approximation to the Willmorc flow with volume 
constraint is established. The existence proof relies on the underlying gradient flow structure of the 
problem: the time discrete approximation is solved by a variational minimization principle. 



1. Introduction 



Let be an open bounded subset of M , 1 < iV < 3, with smooth boundary T. We are interested 
in the following evolution problem 



where the nonlinearity j + cr is a smooth double well-potential (for instance, {j + a){r) = (r^ — 1)^/4), 
u is the outward unit normal vector field to F, and / denotes the spatial mean value of an integrable 
function /, namely, 



As one can easily realize from ([T]) and (E]) by integrating over Q, the mean value of v is conserved 
during the evolution, that is, v{t) = vq. 

The initial-boundary value problem is a phase-field approximation of the Willmore flow 

(cf., in particular, [5l[6]) which belongs to a class of geometric evolutions of hypersurfaces involving 
nonlinear functions of the principal curvatures of the hypersurface. Recall that the Willmore flow 
with volume constraint for a family of (smooth) hypersurfaces (S(t))j>o reads 



(1) 
(2) 
(3) 
(4) 



Vv ■ V 



dtv 



Afi + (j + (t)"{v) /i - (j + aYiv) fi = 0, {t,x)e (0, oo)xn 
-Av + {j + a)'{v) , (t,x) G (0,oo) X fi, 

V/i-z/ = o, (t, x) e (0, cx)) X r , 

vq , X E , 




(5) 
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where V, if, and denote the normal velocity of S, the sum of its principal curvatures (scalar 
mean curvature), the product of its principal curvatures (GauB curvature), and the Laplace-Beltrami 
operator on S, respectively, while A is the Lagrange multiplier accounting for the volume conservation 

/ Vds = Q. 
is 

In addition, the Willmore flow is the L^-gradient flow of the Willmore energy 

(6) £w{^) := I ds. 

Related geometric evolution flows involve more complicated energies such as the Helfrich energy 
and additional constraints, for instance on the area, and are found in the modelling of biological cell 
membranes. We refer, e.g., to |2l El IH El |6l |9] and the references therein for a more detailed description 
of these flows and their applications. To our knowledge, the energetic phase-field approximation ([1])- 
(jlj) has been introduced in [5] in order to describe the deformation of a vesicle membrane under 
the elastic bending energy, with prescribed bulk volume and surface area, a related model without 
constraints being considered in [7]. Here, we restrict our analysis to the case of only the volume 
constraint, leaving the more complex case of two constraints as in [S] to a subsequent investigation. 
A nice feature of already reported in [5] is that it inherits the gradient flow structure of the 

Willmore flow and it is actually a gradient flow in L^{Q) for the functional 

(7) E{v) := ^ /" [-Av{x) + (j + ay{v{x))f dx , 

a property which is a cornerstone of the forthcoming analysis. The connection between the minimizers 
of the Willmore energy ([6]) and those of a suitably rescaled version of the energy (JTl) of the stationary 
phase-field model has been investigated in [H El Ej , and we refer to [SI El E] for the analysis of the 
relationship between the phase- field approach ([I])-(11D and the Willmore flow, with or without volume 
and surface constraints. However, the well-posedness of the phase-field approximation does not seem 
to have been considered so far, and the aim of this note is to show the well-posedness of under 
suitable assumptions on the data: more precisely, we assume that there is Co > such that 

(8) J G C^(M) is a convex function with j(0) = j'(0) = , 

(9) a e C^R) with a" E L°°(R) , 

(10) J + (7 > and r (j + a)'(r) > -Co , rGM. 

Next, owing to the already mentioned expected time invariance of the spatial mean- value of solu- 
tions to ([I])-(11D, for a G M we define the functional space 

(11) W := {w e H^{n) : Vw-z/ = on V} and its subset Wa := {w e W : w = a} . 
The paper is devoted to the proof of the following existence and uniqueness result. 

Theorem 1. Given a G M and vq G Wa, there is a unique solution v to ([I])-(l4]) satisfying 
V eC{[0,T];L\n))r]L°°{0,T;Wa) and := -Av + {j + ay{v) e L\0,T;W) 
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for allT > 0. In addition, 

(12) t I — >■ E {v(t)) := -||/i(^)||2 is a non-increasing function 



(13) / -A/i(t) + (j + a)"{v{t)) fi{t) - (j + a)"{v) fi{t) dt < 2E{vo) . 



2 

Owing to the above mentioned gradient flow structure, a classical approach to existence is to use 
an implicit time scheme and solve a minimization problem at each step, see, e.g., [1] or [ini Chap. 8]. 
The existence of a minimizer to the corresponding stationary problem is discussed in Section [21 
and Subsection 12.11 also collects some properties of the auxiliary variable The time discretization 
is next implemented in Subsection 12.21 and convergence of the time discrete scheme is proved in 
Subsection 12.31 with the help of monotonicity and compactness properties. Finally, uniqueness is 
shown in Section [3] by a standard contraction argument. 

2. Existence 

2.1. The energy functional. Following [5], we define the functional E oxiW hy 

(14) E{w) ■.= \ I [-Aw(x) + (j + a)'{w{x))f dx . 

Observe that E is well defined for any w E W thanks to the continuous embedding of H'^{Q) in 
L°°{n), dH]), and dl]). Indeed, for w eW,we have w e L°°{fl) and 

|(j + anw)\ < r f{r) dr + |a'(0)| + \\a"\UM < k'(0)| + f sup {/} + \\a"\\S] \w\ . 

Jo \[-\\w\\^,\\w\\ac] J 

Consequently, (j + cr)'(w) e L^(f2) and E is well defined. We gather some properties of E in the 
next lemma. 

Lemma 2. Given a G M, there is Ci{a) > depending only on Q, a, Cq in ( ITOj) . and a such that 

(15) \\w\\H2 + \\j'{w)\\2<Ci{a) (i + /eh) for all w e . 

Proof. Consider w G Wa and put /x := —Aw + (j + (7)'(w). Then /x G L^(f2) with \\fJ.\\l = 2E{w), 
and we infer from (fTOD that 



/ w fi dx = \\Vw\\l + w {j + o-)'{w) dx > \\Vw\\l - Cq \Q\ 
Jn Jn 

Combining the above inequality with the Poincare-Wirtinger inequality 
(16) — w^lh < ^"2 II Vw||2 , 
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we obtain 

||Vtf II2 < Co\^l\ + / w fi dx < Col^l + ||w||2||/x||2 
Jn 

< Co\n\ + y^2E{w) + \\w - ah) < CqM + ^2E{w) + \\Vw\\2) 

< CQ\n\+a\n\^'^ V2E{w) + ^ ||Vu;||^ + C2^ E{w) , 
hence || Vw||2 < C{a) (1 + E{w)). Using again ( fT6l) . we conclude that 
(17) \\w\\m<C{a) (l + /BH) • 

Now, w E W solves —Aw + j'iw) = fJ^ — (t'{w) and, owing to the monotonicity of j', a classical 
monotonicity argument shows that 

||A«7||2+||/HI|2< ||/i-a'H||2. 

It then follows from (j9]) that 

IIA^Ih + II/HII2 < M2 + WmM'/' + l|a"lU \\wh , 



which, together with (JT7I) and ||/i||2 = ^/2E{w), gives ( ITSj) . □ 
Next, given r > and / G we define the functional F^-j on TV by 

(18) F,j{w):=^\\w- f\\l + T E{w), w e W . 

Lemma 3. Given a G M, the functional Frj has (at least) a minimizer in Wa- 

Proof. We set F := F^-j to simplify notations. Since E is nonnegative, F is obviously nonnegative 
and there is a minimizing sequence {wn)n>i in Wa such that 

(19) m« := inf {F(w)} < F(w„) < m„ + - , n>l. 

Since F{wn) > t E{wn), we readily infer from ( 1T^ that (i?(w„))„>i is bounded, a property which 
in turn implies that {wn)n>i is bounded in H'^{Q) by Lemma [2J Owing to the compactness of the 
embedding of H'^iVt) in C(f2), we deduce that there are w G H^{Q) and a subsequence of (wn)n.>i 
(not relabeled) such that 

(20) Wn — > w in C{Q) and ty„ ^ in H'^{Q) . 

Clearly, the first convergence implies that ((j + o')'(w„,))„>]^ converges towards (j + ct)'(w) in L^(f2) 
and therefore 

F{w) < liminf F{wn) < rria ■ 
As w obviously belongs to Wa by fl20|) . we also have F{w) > rria and is a minimizer of F in VFq,. □ 
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We next derive an energy inequality and the Euler-Lagrange equation satisfied by minimizers of 
Frj in Wa when f = a. 

Lemma 4. Consider a G M and a minimizer w of Frj in Wa- Assume further that f = a. Then 
/i := —Aw + (j + belongs to W , 

'w- f 



T 



-Afi + (j + a)"{w) fi-{j + a)"{w) fi 



ip dx = for all tp e W 
w-fh 



< 



(21) / A/i+(j + a)"H + /i 

and 
(22) 
Proof. We set 

/i := -Aw + (j + (t)'(m;) . 

Consider e G (0, 1) and (p G Wq. As w + e(p belongs to Wa, we have F^-j^w) < Ft-j{w + Sip) from 
which we deduce by classical arguments (after passing to the limit as e — )■ 0) that 

- I {w - f) ip dx + I fi {-Aip + (j + cr)"(w) ip) dx>0. 

Since the above inequality is valid for ip and —(p, we actually have the identity 



(23) 



r 



— {w — f) (p dx + jj, {—Aip + {j + a)"{w) ip) dx = 



for all ip G Wq. Now, if ^/^ G W, the function tp — tp belongs to Wq and it follows from ( l23l) that 
1 



(24) 



-r Jn 



{w - f) dx + / fi {-Atp + (j + o-)"{w) tp) dx = (j + (r)"{w) fi dx 



since w and / have the same mean value a. Since fi G L^(fi) solves the variational equality flMj) for 
all test functions G VT, we deduce that fi E W and satisfies ( 12T|) . 
Next, for r/ G (0, 1), let (yj^ be the unique solution in Wq to 



iPr,-ri Aiprj = -Aji + (j + 0-) (u;) /i - ( j + cr)"(u;) /i in , 

the right-hand side of the previous equation being in L^(f2) since fi E W and w G H'^{Q) is bounded. 
Also, the right-hand side of the previous equation has a zero mean-value so that iprj G Wq- Taking 
ip = ipr] in f l2T]) . we realize that 

- / 



+ iprj-r] Aipn 



iPn dx = Q . 



from which we deduce that 



\Wn\\l<\M\l + V llVy^.l 



w-/ , ^ Ik -/lb II II 



T 
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whence 

II II ^ Ww-fh 



Since {(Prf)r] converges toward (— AyU+ (j + cr)"(ty) /i — (j + a)"{w) //) in L^(f2) as r/ — )■ 0, f l22|) follows 
from the above inequality. □ 

2.2. Time discretization. Let a G M and take an initial condition vq G W^. We consider a positive 
time step r G (0, 1) and define a sequence (f^)n>i inductively as follows: 

(25) v^,:=Vo, 

(26) f^^i is a minimizer of -Ft,d^ in VF^ , n > , 
the functional -Fr,i);; being defined in ( ITH]) . Setting 



(27) ;,;:=_A< + (j + a)'«) and M^= (j + a)>;;) /i^ , 
we define three piecewise constant time-dependent functions v'^, fi"^, and M'^ by 

(28) {v^{t), /i"(t), M"(t)) := «, /i^, M;) for t G [nr, (n + l)r) and n > . 
Lemma 5. For r G (0, 1), ti > 0, and t2 > ti, we have 

(29) E{v^{t2))<E{v^{h))<E{vo), 

(30) \\v^{h)-v^it,)\\l<2E{vo) (r + h-h), 

||-A/i-(t) + (j + a)"(t;^(t)) /i^(t) - M-(t)||^ rft < 2E(t;o) . 
Proof. Consider n > 0. Since G PVq,, we infer from fl26|) that F,- j,r (t;^_^-^) < F^-^yr^v^), that is, 
(32) ^ + <i?«) . 

Let ^2 > ^1 > and put := [tj/r] (the integer part of tj/r), z = 1,2. On the one hand, 71-2 > 'n-i 
and it readily follows from ( 132|) by induction that 

E(.;-(t2)) = i?KJ < i?KJ=i?(i;-(ti)), 
whence (l29l) . In particular, we have 



(33) ^ sup ||/i^(t)||^ = snpEiv^it)) = supE «) < E (t;J) = E(t;o) . 

^ t>0 t>0 n>0 

On the other hand, summing (132|) over n G N gives 

^ oo 

(34) - IK+i - <\\l <E(<) = EM , 
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from which we deduce that 



n2 — 1 



n=ni 



1/2 / \- |L,r _ r||2 
+1 ^n\\2 



< + i2rEivo)f' 
and thus fl30l) . Finally, for > 0, we have v^_^-^ = = a by fl26l) and we infer from fl22|) that 

K+l - <l|2 



\n=ni 
1/2 



+ (j + - ^:+i||2 < 



r 



Combining ( !34|) and the previous inequality give 



oo 

||-A/i^(t) + (j + /i^(t) - M-(t)||^ 

°° r{n+2)T 

< E / ll-^/^n+i + (j + z^^+i - dt 

n=0"^("+l)^ 
°° II r '''112 

< J2 "'^^^ <2Eivo), 



T 
n=0 

and the proof is complete. □ 

Useful bounds on {v'^)r and (/i^)T follow from Lemma [51 
Corollary 6. For all T > 0, there is C^{T) > depending only on a, Vq, j , a, and T such that , for 

r G (o,i)n(o,r), 

(35) sup Wv^m^, < C3(T), 

te[o,T] 

(36) £ {Wl^'mm + Wf^^mU dt < Cs{T). 

Proof. The boundedness (!35l) of (f^)r is a straightforward consequence of (ITS!) and (!33|) . Next, owing 
to the continuous embedding of H'^{Q) in L°°(f2) and fl35l) . the family ((j + (T)"(f'^))T- is bounded in 
L°°((0,T) X fi) which, together with (jSS]), imply that 

(37) {{j + a)"{v^)fi^)r is bounded in L'^{0,T; L\n)) . 
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Setting := -A/i^ + (j + ct)"{v^) fi^ - M\ it follows from (EI]) and ([37]) that 

l|A/i^(t)||^ dt^ ' = (^J^ WU + ariv^it)) ^^^{t) - M-{t) - r{t)\\l dt ^ ' 

T \^/^ / rT \ 1/2 

2 J , \ , / / II xr/,M|2 



< 2 i^j \\{j+anv^{t))f,^mtdtj wrmtdt 

which gives the boundedness of (yu"^),- in L2(r,T; H^{n)) with the help of ([33]). Finally, fi^ e W and 
solves 

-A/i" + j"{v^) = r - a"{v^) fi^ + M^ in Q. 

Taking the scalar product in L^(f2) of the previous equation with /i"^ and using the nonnegativity of 
j" due to the convexity ([H]) of j and the boundedness ([9]) of a", we obtain 

l|V/i1|^ < + [ f{in il^n'dx 

Jn 

< wrh ii/iii2 + ik"iioo \\i2^\\i + \M^\ wi^^h- 

We next deduce from (15^ and (1571) that 

wvfi^wi <c{T) (i + iirib) , 

and the boundedness of the right-hand side of the above inequality in L'^{t,T) follows at once from 
([HI]). □ 

2.3. Convergence. Owing to ( 130]) . ( l35l) . and the compactness of the embedding of H'^{Q) in C{Q), 
a refined version of the Ascoli-Arzela theorem (in the spirit of P, Prop. 3.3.1]) ensures that (f^)r is 
relatively compact in C([0,T] x Q) for all T > 0. Consequently, there are three functions v, fi, and 
M and a subsequence of (f'^)r such that, for all T > 0, 

V eC{[0,T]xn)nL^{0,T;H\n)) , ^ G L°°(0,T;L2((])) , MeL°°(0,T), 

and 

(38) v^"^ — > V in C([0,r] x H) , 

(39) v^^ ^ V in L'^{0,T;H\Q)) , 

(40) ^ /i in L°°(0,T;L2(fi)), 

(41) M^^ ^ M in L°°(0,T). 

Thanks to the smoothness of j and a and the convergences ( l38l) - (l^Il) . it is straightforward to pass 
to the limit in (1271) and conclude that 



(42) /i = -Av + (j + cr)'(t;) and M = ( j + cr)"(tO . 
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In addition, ( l36l) . ( I40l) . and a lower semicontinuity argument guarantee that 

(43) jie L\0,T;H\n))nL\0,T;H\n)) for all T>0. 

It remains to derive the equation solved by v. Let ip G W, t > 0, n = [t/r], and m G {0, . . . , n — 1}. 
Using the definition of v'^,-^ and Lemma HI we are led to 



n 



-71 ' ' 

^rn at i / ■ i r \ r 71 /rr 

A/i^+i + (j + c^) f^m+i - 



ip dx = 



which also reads 

/■(m+2)r /• 

K+1 -vl) i,dx= / [A/x-(s) - (j + a)"{v^{s)) fi^is) + M-(s)] ^^J dxds . 

'n J{m+i)T Jn 

Summing the above identities over m G {0, . . . , n — 1}, we obtain 

l>(n+l)T n 

{< -vl) ^dx= / [A/i"(s) - (j + (y)"{v^{s)) + M"(s)] c^xrfs . 

Jt Jn 

/■(n+l)r 

{v^{t) -vo) ip dx= / [A/i"(s) - (j + fT)"(t;"(s)) + M"(s)] ^ rfxds . 

Jt Jn 

Noticing that t < {n + l)r < t + r, we may take r = in the above identity and pass to the limit 
as /c — )■ 00 with the help of (!38!) - (l4T!) to obtain 

(44) [ {v{t)-vo) %l)dx= f [ [Afi{s) - {j + a)"{v{s)) fi{s) + M{s)]iIj dxds . 

Jn Jo Jn 

Collecting ( I42l) -( l44|) completes the proof of the existence part of Theorem [H The properties (fT2l) and 
(fT3|) next follow from ([21]), (jSl]), and the convergences (l^-(|4T|). 

3. Uniqueness 



Let Vi and V2 be two solutions to ([I])-(ll]) with /ij := —Avi + (j + o")'(f j) and Mj := {j + a)"{vi)fii, 
i = 1,2. Fix T > 0. Since H^{Q) is continuously embedded in L°°{Q), the regularity properties of 
^1, V2, fii, and /i2 listed in Theorem [1] ensures that there is K > depending on T such that 

(45) sup (||t;i(t)||oo+IK'2(t)||oo + ||/ii(t)||2 + ||/i2(t)||2)+ / ( || /ii (^^) || L + II (s) || L) ds < K . 

te[o,T] Jo 

It then follows from fH5l) and the smoothness of j and a that 

(46) + + 

< IKj + ^)"'IIl-(-av<) 1^1-^21 l/^il + IKj + c^)"IIl-(-k,x) I/^i-/^2|, 

< C (|/ii| - t'sl + l/ii - Ai2|) , 
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from which we deduce that 

(47) IM1-M2I < jjU + (^n^,) fi, - {j + aYiv^) fi2\ dx 



< C \vi - V2\ + IfJ^i - fi2\) dx 

Jn 

< C (||/il||2 - ?;2||2 + Wl^l - /^2||2) • 

Since Wi — V2 solves 

dt{vi - V2) - A(/ii - /is) = Ml - M2 - (j + (t)"{vi) /ii + (j + (J)"{v2) 122 
and Vi — V2 and /ii — fj,2 both belong to W, we have 

^ ^11^1 ~ ^'2 II 2 = / - 1^2) ^{vi - V2) dx+ [ (Ml - M2) {vi - V2) dx 
2 dt Jn Jn 

[U + (^y'i^i) - (j + ^)"(^'2) /i2] (f 1 - f2) dx . 
We deduce from and ^ that 

2 ^IK'l - ^'2||2 = / (/^l - /^2) [(j + f^)'(^l) - (j + C^)'(f2) - (/il - /i2)] 



+ / (Ml - M2) {vi - 1^2) c?a; 

[(j + o-)"(t'i) /ii - (j + o-)"(t;2) /i2] (t'l - V2) dx 



< 


IKj + f^)"llL-(- 




il - /i2||2 ll^^l - ^^2||2 


- 11/^1 -/^2||2 


+ 


C (ll/iilh IK'i - 


- ^'2||2 ^ 


- 11^^1-/^2112) ll^^l - 


^^2||2 


+ 


c / (I/ill K'l 


-t'2| + 


l/Wl - At2|) \Vl - V2\ 


dx 


< 


C ll/il - /i2||2 1 


\Vl - V2 


II2 - Wf^l - f^2\\l 




+ 


C (1 + ll/illloo; 


) hi - 


V2\\l 




< 


c (i + ii/xiiu; 


) \\vi - 


V2\\l. 





Therefore, recalling ( H5|) . 

11(^1 - ^2)(t)||^ < ||(t;i - t;2)(0)||^ exp (1 + ||/ii(s)|U) rfs ) < C ||(^;i - ^;2)(0)||^ 

for t G [0,T], and the uniqueness assertion follows. 
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